Stochastic resonance and the trade arrival rate of stocks 
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Abstract 

We studied non-dynamical stochastic resonance for 
the number of trades in the stock market. The trade 
arrival rate presents a deterministic pattern that can 
be modeled by a cosine function perturbed by noise. 
Due to the nonlinear relationship between the rate 
and the observed number of trades, the noise can 
either enhance or suppress the detection of the deter- 
ministic pattern. By finding the parameters of our 
model with intra-day data, we describe the trading 
environment and illustrate the presence of SR in the 
trade arrival rate of stocks in the U.S. market. 



1 Introduction 

It can be very useful to model and understand mar- 
ket microstructure. For instance, it has been shown 
that the number of trades and volume can help ex- 
plain volatility fluctuations and therefore help model 
asset returns [^ (TUl [55]. More importantly, these 
studies link market activity (and therefore market 
microstructure) to the large class of asset allocation 
and pricing models based on time changed Brownian 
motion [9l [23l [27| . More practically, we also point 
to the increase of intra-day trading that some claim 



should account for more than 50% of the trading vol- 
ume by 2010. High-frequency trading will increas- 
ingly require market microstructure for risk control 
and trading. 

In this paper, we use the doubly stochastic Pois- 
son process to model the number of trades inside of a 
trading day. Using this model, we focus on the inter- 
action between noise and the trade rate "U" -shaped 
pattern (higher average arrival rate at the open and 
close of the trading day). We report an interesting 
feature of such model: stochastic resonance (SR). 
The presence of stochastic resonance indicates that 
a correct amount of noise is, in fact, beneficial to in- 
crease the detection of the "U" -shaped pattern. That 
is, random trading (represented by random arrival 
rates) can constructively contribute to the existence 
of such periodic patterns. 

1.1 "U" -shaped pattern 

Finance and economics show a variety of well doc- 
umented periodic patterns [T3]. In particular, there 
is the intra-day "U" -shaped pattern (or patterns). |1| 
That is, the average intra-day volatility, the trad- 



^Also sometimes called the intra-day "smile" or inverse J 
pattern. 



ing volume, bid-ask spreads, trading activity (num- 
ber of trades), and even returns are significantly 
higher near the open and close of the trading day 
[3 [201 [23 m [3D]. These quantities form on average 
an "U" shape; hence, the name " 'U'-shaped pattern" . 
In fact, it has been shown that there need not be 
a official exchange halt. Studies using 24-hour mar- 
kets (such as currency markets) find that "U" -shaped 
patterns are also present within the normal business 
hours PI]. 

There is a large literature that explains the dif- 
ferent "U" -shaped patterns. In general the litera- 
ture has been grouped into two lines: The first line 
suggests that asymmetric information is a necessary 
condition for the "U" shape to exist. By asymmet- 
ric information, it is meant that some traders have 
private information about the future asset returns [T] 
or payoffs [18]. It is the interplay between informed 
traders and uniformed traders together with market 
stops that will lead to the observed pattern both in 
returns and volumes [TB]. 

The second line contests the importance of asym- 
metric information by pointing out to the presence 
of the same patterns in markets with no asymmetric 
information [TTl [12] . They argue that the presence of 
market stops could be enough. This line emphasizes 
the exchange of risk close to market stops [TB] . The 
argument is that information accumulated overnight 
will generate an increase of the trading activity at 
market open because portfolio managers will want 
to bring their portfolio to a particular risk target as 
soon as possible. The same will happen at the end 
of a trading day, when traders will want to prepare 
themselves for the next overnight period. Therefore, 
the pattern could be just due to the timing of trans- 
action decisions by institutional investors and long- 
term investors with short term investors representing 
the other side of the trade [5] . 

In this paper we do not address the micro- 
economical origin of such patterns as, for instance, 
in Hong and Wang [TH]. We take an empirical ap- 
proach to the problem. We start with a plausible 
statistical model which describes the data well. Since 
our model presents stochastic resonance, we conjec- 
ture that such a phenomenon can be important to 
the "U" pattern. The presence of SR could indicate 



that noise traders can, in fact, magnify such pattern. 
However, in the present paper we will not be able to 
reach such level of description, leaving such questions 
for future studies. 

From a practical point of view, our work still 
presents useful results. For instance, one can also 
think on a possible algorithm where random trading 
(hence undetected to other traders) can enhance even 
further such a pattern. Of course, the benefit of en- 
hancing such a pattern is, a priori, not clear. But in 
a similar fashion there could be other patterns where 
similar analysis will lead to undetected (and bene- 
ficial to the trader) market manipulation. We also 
point to the analogy drawn recently by Dempster et 
al. [13] as a possible starting point to generalize SR 
to asset allocation, for instance. 



1.2 Stochastic resonance 

Stochastic resonance (SR) has been successfully used 
to explain signal transmission and detection in a wide 
variety of systems in diverse areas such as physics, ge- 
ology, neural science, engineering, and biology [T5] . 
SR is briefly defined as noise-enhanced signal detec- 
tion. Typically there is a threshold below which the 
signal cannot be detected. Adding noise to the sig- 
nal can increase the number of events that cross such 
a threshold, hence increasing the detection of such 
signal. Imagine a sine wave of amplitude less than 
the detection barrier. Noise added to the sine wave 
will increase the events that cross the detection bar- 
rier. However, if the noise intensity is too high, it will 
corrupt the signal beyond recognition. 

In finance, stochastic resonance has been suggested 
to explain market crash and high-frequency patterns 
in the exchange rate [TU [25 . Different from previous 
work, our paper employs non-dynamical, threshold- 
free stochastic resonance [4j [5]. Non-dynamical SR 
is purely stochastic in nature and, therefore, it does 
not require a microscopic model (such as Newton's 
equations of motion) . Consequently, it can be more 
general because it dissociates a specific microscopic 
model from the SR phenomenon. To the best of our 
knowledge, this is the first time this particular SR 
class has been applied to financial data. 



Non-dynamical 
resonance 



stochastic 



In direct analogy to non-dynamical SR in other dis- 
ciplines, we assume that trade arrivals follow a time- 
dependent Poisson process. We take the point of view 
of a market participant and observer. We observe 
signals in the market that are being transmitted to 
all market participants. Our goal is to describe the 
present state of the trading environment and there- 
fore detect the presence of SR. 

Non-dynamical and threshold-free stochastic reso- 
nance was originally discovered for systems that can 
be modeled with a doubly stochastic Poisson process 
(or time-dependent Poisson process) 5 . Such a pro- 
cess is constructed out of Poisson process N{t) with 
the probability distribution of n events during time 
interval T given by 



-rT 



PN{n) = 



{rTY 



(1) 



where the event arrival rate r is itself a time de- 
pendent random process independent of N{t) [3 . 
Bezrukov and Vodyanoy have shown that the follow- 
ing arrival rate r(i) of a doubly stochastic Poisson 
process exhibits SR: 



r(i)=r(0)exp(yl(t)), 



with A{t) as 



A{t) = W{t) + As cos{2tt fst). 



(2) 



(3) 



That is, A{t) is the sum of zero mean Gaussian noise 
W(t) with root mean square a; and a weak periodic 
signal with amplitude A^ and frequency fs [5]. 

In the case of Equation ([3]), the slow cosine- wave 
As cos(27r/si) can be better detected, due to SR, by 
analyzing the final Poisson generated time series N{t) 
if there is a correctly tuned noise W{t). However, 
excessive noise might corrupt the signal. 

2.1 Signal-to-noise ratio 

For a quantitative analysis of the synchronized pro- 
cess, signal-to-noise ratio (SNR) is introduced to 



measure stochastic resonance. SNR is defined as the 
height of the signal divided by the height of the noise 
background of the power spectrum at the signal fre- 
quency fs [13] . We speak of SR if there is a maximum 
in the output SNR as a function of noise intensity a. 
That is, there is an optimal value of the noise that 
concurs with the periodic signal. 

Bezrukov and Vodyanoy derived the output SNR 
for the doubly stochastic Poisson process with the 
rate given by Equation ^ and Equation ([3]) via 
power spectral density for As <^ a and fs much 
smaller than all other characteristic frequencies |3 
The output SNR is given by 



SNR 



(A2r(0)/2)exp(aV2) 



2 + (2r(0)A/,) exp(aV2) E^=i(l/")(<^'"/"!) ' 

(4) 
where fc is the noise corner frequency [H \5\ . We also 
have fc = (27rrc)~^, where t^ is the noise W{t) corre- 
lation time. The noise has a Lorentzian power spec- 
trum with exponential autocorrelation function. We 
have also verified Equation ^ by Monte Carlo (MC) 
simulating Equation ([1]) with rate given by Equa- 
tions ([2|) and ([3|) 8 . Using our MC simulations, we 
remark that Equation (j4]) is still a very good approx- 
imation outside the range of parameters for which it 
has been derived [S]. 

The maximum SNR as a function of a in Equa- 
tion (|4]) is the result of a competition between the 
exponential in the numerator and the exponential 
times a linear term (to first order approximation) in 
the denominator. The constant j4^r(0)/2 is a general 
scale of the amplitude of the maximum. The loca- 
tion of the maximum is given by r(0)//c. r(0)//c 
is the number of events that occur over the correla- 
tion time. If r(0)//c is large, the maximum moves 
to small a; if r(0)//c is small the maximum moves 
to large a. Therefore, the existence of the maximum 
depends entirely on r(0)/fc: SR is more important 
for a random variable where the statistics is poor; 
that is, the equilibrium arrival rate r(0) is small. SR 



^They have, in fact, derived it for a sine-wave signal; how- 
eve, the power spectral densities for sine and cosine are iden- 
tical. Using cosine instead of sine here will not change the 
SNR. 



is also more important when the noise present (or 
added) to the system is uncorrelated. Furthermore, 
it also shows that no matter how good r(0) is, the 
output signal detection could be further improved if 
fc is big enough. Equation ^ shows that if the con- 
dition fc > 2r(0)7r [Tcr(O) < 0.25] holds, the process 
demonstrates SR. 

3 Numerical experiments 

We cannot always perform controlled experiments 
in economics; therefore, we need to extract all the 
needed information from the data we observe. We 
first parameterize the time dependent rate r{t) de- 
fined by Equations ^ and ^ without the Gaussian 
noise W{t) by assuming that parameters for the rate 
are constant for at least 1 month. We have 



r(t)=r(0)e^^'=°''(2^^»*\ 



(5) 



where fs is one over the trading day, fs — 1/23400 
seconds. Notice that once we have detected the pe- 
riodic signal As cos{2tt fst), we have bypassed one of 
the potential benefits of SR: increasing statistics to 
help the signal detection. Therefore, we cannot infer 
if our signal reconstruction was helped by SR. How- 
ever, we can still say that a system can potentially 
present SR by identifying all the parameters. 

In our numerical experiments, we obtained tick- 
by-tick financial data (smallest time resolution is 1 
second) for 12 trading stocks (AA,AEP, AES, AMG, 
AOL, ATI, EK, F, JPM, PFE, PG, and S) during 21 
trading days (January 2001). First, we analyze the 
daily pattern by constructing a proxy for the trade ar- 
rival rate inside of the 6.5 hours composing the trad- 
ing day. We bin the data in 2-minute (120-second) 
intervals using 1-second resolution from the first day 
to the last day (total of 21 days). We assume that 
the rate in each of the 2-minute bin is constant. We 
remove the first 5 minutes at the beginning of the 
trading day due to data inaccuracy. We end up with 
a total of 192 2-minute intervals in a day. 

We then calculate the arrival rate for every 2 min- 
utes by counting the number of trades in every inter- 
val and dividing by 120 seconds; that is, the average 
number of trades in each of the 192 bins during 21 
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Figure 1: Normalized average intra-day arrival rate pattern for 
12 stocl<s. The number of trades per unit of time (2-minute bins 
are used here) are chosen as a good proxy for the arrival rate of 
trades. The symbols are the average over 21 days and 12 stocks. 
The error bars represent 1 standard deviation of the empirical 
data. The solid line is the theoretical model in Equation Js}. 



days. This is equivalent to fitting a constant rate 
Poisson process since the maximum likelihood esti- 
mate for the rate is simply the average number of 
events per unit of time. Once we have the empirical 
daily rate pattern, we fit the daily rate pattern to 
Equation ^ to estimate r{0) and Ag. 

FIG. [T] shows the arrival rate averaged over 21 days 
and all 12 stocks. In order to average over all stocks, 
we subtract ln(r(0)) from the natural logarithm of 
the empirical rate of each stock and then multiple 
by 1/As to obtain cos{2tt fst). Thus, we are able to 
recover the parameter independent exp(cos(27r/st)). 
Equation ((51) is a very simple periodic function and 
hence the most simple example one can choose. How- 
ever, FIG.[l]shows that it gives a good approximation 
to the empirical daily rate pattern especially for the 
second half of the day. The first half is on average sys- 
tematically higher indicating that a skewed function 
could be better. Nonetheless, since r(0) measures 
the vertical shift of the whole curve, the value of r(0) 
will not change very much even if asymmetry is taken 
into account. The quality of our approximation will 
suffice for defining SR. 

After r(0) and As have been estimated, the remain- 
ing fluctuation are assumed to be generated from the 
random process W{t) with intensity a [Eq. ^]. At 
this stage, we only need to find Tc (or fc) [Eq. (j4|)] to 




Figure 2: Autocorrelation of the number of trades for 2 stocks 
(ATI and F). The symbols represent the empirical autocorrelation 
for the data, and the solid lines represent the theoretical autocor- 
relation constructed via MC simulation. The inset presents the 
behavior of the tails in log-linear scale. 
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Figure 3: The top panel of the graph shows the daily trading 
rate pattern [Eq. ((S}] for F. The bottom graphs are the selected 
PDFs of 2-minute intervals along the trading day for F; the sym- 
bols are data points, and the lines are the theoretical PDF. Once 
the daily pattern is found (top panel), we fit all 2-minute empirical 
PDFs along the day with one single parameter a (lower panel). 



demonstrate the presence of SR for a given stock. 

We find Tc indirectly through the autocorrelation 
function of the number of trades iV(t). Notice that t^ 
measures the characteristic memory of the Gaussian 
noise W(t) in Equation ^. However, we cannot mea- 
sure the autocorrelation of W{£) since Wif) cannot 
be observed from the data directly. Furthermore, we 
can achieve a 1-second resolution by using the auto- 
correlation N(t), since iV(i) is directly observed from 
the market. 

We find the autocorrelation of the number of trades 
for our model by Monte Carlo (MC) simulating Equa- 
tions H)), (l2|), and ([3]). We simulate a time series of 
the number of trades N{t") with a band-limited white 
noise W{€) of corner frequency /c. For every /c (or 
T^ there will be a time series of the number of trades. 
For every such time series, we have an autocorrela- 
tion function. The correlation time t^ for a given 
stock is the one which minimizes the mean square 
error between the simulated autocorrelation and the 
empirical autocorrelation function. 

FIG. [D presents the autocorrelation of the number 
of trades Nify for two stocks, ATI and F. We calcu- 
late the autocorrelation function using all 21 days of 
tick- by-tick data (a total of « 530000 data points). 
These stocks represent the typical autocorrelation for 
our 12 stocks. We are able to model autocorrelation 
function up to correlations of the order of 1%. How- 
ever, in agreement with the literature, some stocks 
clearly present a slow decaying autocorrelation. No- 
tice that F has a slow decaying tail above the au- 
tocorrelation noise level (« 3/v^530000) = 0.0041). 
ATI, however, decays faster and has no significant 
autocorrelation. 

The presence of a persistent tail for F shows that 
the autocorrelation could better modeled either by 2 
exponentials (two different Tc) or by a power law. The 
exponential autocorrelation function with a single t^ 
is, therefore, a first order approximation. However, 
such first order approximation will suffice to define 
SR. We note that is highly unlikely that t^ will be 
small enough for a stock to present a maximum in 
SNR and still have a large and significant autocor- 
relation tail. That is, most of the autocorrelation 
(> 99%) is generally explained by a single exponen- 
tial function. If a stock has a large autocorrelation 



tail, it will have a large Tc as happens with F (and 
not with ATI). 

The noise intensity a defines the present trading 
environment in the SNR graph [FIG. (g])]. We find 
(T by fitting the theoretical probability density func- 
tion (PDF) of the doubly stochastic Poisson pro- 
cess to the empirical PDF j^S]- Note that r(0) and 
As were estimated previously via Equation ([5]). We 
replace the previously estimated r(0) and As into 
Equation ^. We continue with the approximation 
that every 2-minute interval is described by doubly 
stochastic Poisson process with constant arrival rate 
r(0)exp(AsCOs(27r/si)) perturbed by noise. Thus, 
within each 2-minute interval, the noise can be de- 
scribed by the log- normal distribution exp(iy(i)). 
We apply numerical integration to calculate the in- 
tegral of the log-normal (which defines the random 
rate) and the Poisson (which takes this random rate) 
to find the PDF of number of arrivals in every 2- 
minute interval. Since the number of arrivals in 2- 
minute intervals are small, we compensate the small 
number of arrivals by fitting all the 2-minute inter- 
vals in a trading day at once to obtain a by the mean 
square method. 

FIG.[3]shows the daily trading pattern by binning a 
trading day into 2-minute intervals for F (other stocks 
present similar results) . The top part of the figure is 
the daily trading pattern. For a cosine function, the 
first half of the day is the mirror image of the second 
half. Wc put both halves on top of each other. The 
solid thick line is obtained from our model. After 
finding the parameterized mean, all the data points 
are fitted at once to estimate the noise volatility a. 
The bottom part of the figure shows selected PDFs 
for few 2-min intervals along the trading day. For 
instance, the first PDF represents the first and the 
last 2-minute intervals of the trading day. FIG. [3] 
shows that our simple model with a single a can fit 
the empirical PDF up to 4 orders of magnitude. This 
remarkable agreement is found for all stocks studied 
along the entire trading day. 

Once a is obtained by our scheme, we have de- 
scribed the trading environment using our model; 
thus, stocks that present SR can be identified. FIG. 
|4] shows the SNR as a function of noise intensity. 
We have normalized all SNR curves such that the 
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Figure 4: Output signal to noise ratio (SNR) of 4 trading stocks 
(AEP, ATI, PG, F) using Eq. JJ]). The circles are the current value 
of cr. The figure shows that SR occurs at different levels in AEP 
and AIT, but not in F and PG. 





r(0)(l/sec) 


As 


CF 


Tc{sec) 


Tcr(O) 


AEP 


0.0283 


0.2880 


0.9554 


1.7000 


0.05 


ATI 


0.0146 


0.1253 


0.6688 


4.7000 


0.05 


PG 


0.0946 


0.2856 


0.6839 


2.6000 


0.26 


F 


0.0947 


0.2194 


0.8118 


26.5000 


2.60 



Table 1: Parameter values for FIG. 4. 

SNR for cr = is equal to 1 for all stocks. Table 
[1] has the recovered values for the parameters of our 
model. We can see that, in fact, stocks that satisfy 
r(0)rc < 0.25 such as AEP and ATI present a max- 
imum on the SNR (and, hence, show SR). From our 
sample of 12 stocks, 1/2 show SR. The largest r(0) for 
the stocks that present SR is w 4 trades per minute, 
whereas, theoretically, we need to have r(0) < 15 
per minute (assuming that Tc = 1 sec is the lowest 
possible value). These results indicate that the less- 
traded stocks are more likely to present SR. Usually 
the less-traded stocks are stocks with smaller market 
capitalization. Therefore, SR should be more impor- 
tant for small cap stocks. 

The values in Table[T]together with FIG.[3]describe 
the statistical trading environment. From the point 
of view of a market observer, the signal of interest 
[As cos{2tt fct)) is being transmitted more efficiently 
(ATI and AEP) or less efficiently (F and PG) with 



the help of noise present in the trading environment. 
From a point of view of a market participant, both 
AEP and ATI are not at the SNR maximum. There- 
fore, adding noise by adding trades with a random 
arrival rate should increase the detection of the ar- 
rival rate "U" -shaped pattern. In the case of AEP, 
for instance, a 25% increase in the noise intensity is 
beneficial to the periodic pattern. 



4 Conclusion 

We have proposed to model the number of stock 
trades inside of the trading day with a doubly 
stochastic Poisson process. With a doubly stochastic 
Poisson process, we are able to model both the deter- 
ministic periodic trend ("U" -shaped pattern) in the 
trade arrival rate as well as deviations from it due 
to noise. We show that such model describes well 
the empirical data and that non-dynamical stochas- 
tic resonance can be a relevant phenomena for the 
"U"-shaped pattern detection. SR as presented in 
this paper is defined as a weak regular oscillation that 
may be amplified by adding noise to the input signal. 
That means that the increase of noise trading (rep- 
resented here with an increase in the random trade 
rate) can, in fact, augment the relevance of the exit- 
ing "U" -shaped pattern. 

We envision that extensions of the model presented 
in this work could be possible for different economical 
time series. In our theoretical framework, the pres- 
ence of a periodic pattern is a necessary condition 
for SR to be relevant. We know that such periodic 
patterns are, in fact, present in economics in general 
(see [m [171 [221 [29] for additional examples), there- 
fore making SR a potentially important phenomena. 
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